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Abstract—A numerical model on the flexibility method in the case of a multilayer beam finite element has been developed 
and the contributions to its recent developments being made at Mechanical laboratory, Department of physics, Faculty of 
Sciences Rabat (Morocco). The results of the experiments and those of numerical calculations were concordant in the 
case of quasi-static loading. These results were based on the approach "finite element" coupled with a non-linear model 
[23]. Firstly, we present here the results based approach "finite element" related to the analysis of a bending square plate 
under concentrated and uniform load, clamped or simply supported on the contour. On the other hand, we present some 
results which we evidence to the problem related to the shear locking. The numerical model is based on a three- 
dimensional model of the structure seen here as a set of finite elements for multilayered plates multi cellular matrix 
(concrete) and a set of finite element fibers for reinforcement. The results obtained confirm the ability of these tools to 
correctly represent the behavior of quasi-statics of such a complex system and presage the deepening of a digital 
tool developed. 
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I. INTRODUCTION 

The phenomenon related to blockade by shear (or appearance of a parasitic stiffness) is a numerical problem that 
drew attention of many researchers in the past twenty years and an abundance of solutions which has been discussed in [3, 9, 
10, 11, 12, 19. 20, 22], One way to avoid the appearance of shear locking and thus make the solution independent of the 
slenderness ratio (the ratio of length L I thickness h) is to calculate the terms oi ihc stiffness matrix by integrating accurately 
the relative terms bending and sub-integrating the terms relating to shear [4,5,6,8,13,14,15,16,17,21 ,22], To improve this 
phenomenon related to the numeric;*! computation and propose a more efficient solution, we developed a model based on the 
flexibility method [23]. The model is formulated on the basis of the forces method by an exact interpolation stresses [18]. 
This makes it possible to calculate the flexibility matrix, which is the inverse of the stiffness matrix. The purpose of this 
study is the modeling of the structural response of the sails carriers subjected to seismic effects using a comprehensive three- 
dimensional numerical model using a nonlinear finite element approach coupled with a damage model developed for the 
behavior of con erete material. In this second paper, drawing on the results of the first article unci those of [1,2 ,7], we present 
only some results related to the analysis of a homogeneous square plate in bending subjected to a concentrated and uniform 

II. MODELING 

Complementary to the trials and their interpretation, numerical modeling of this situation type has several 
advantages. In this case, it already developed an ambitious and effective model capable of taking into account the different 
aspects of this complicated problem, including the quasi static and dynamic loading. Then after this satisfactory model, it 
has to constitute a way to complement the experimental measurements by providing new data. As such, it should contribute 
to a better understanding of the phenomena involved and to further provide a basis for dimensionality development methods. 

1. METHODOLOGY 

An immediate challenge before addressing the simulation of such problems is to choose the right methodology. 
The philosophy retained here is to realize the contribution of research in civil engineering to respond in a context of 
operational engineering. The choice was made on the use of finite element plate's multilayer multistage three nodes and two 
degrees of freedom per node. 

A realistic numerical prediction of the structural response of such a structure requires a rigorous three-dimensional 
geometric model of the system components. This model and its numerical analysis are implemented in the finite element 
code RE-FLEX. 

Then, the plate is meshed by including its geometry in a full mesh adapted to the different areas of the problem (it 
is discredited into layers and its thickness h in cells along x and y the surface) [Fig.l], 
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Figure 1 - Finite Element Model: Efforts resulting in a plate 



Where N xx ,N yy represent the normal forces and N xy the shear plane. M xx , M yy represent the bending 
M torque. T x , T are the transverse shear stresses. 



2. CALCULATING THE ELEMENTARY FLEXIBILITY 

The exact interpolation functions are obtained by writing the various external forces of any point of the finite 
element, which here are the internal forces of the structure, according to the nodal reduced effort. Thus, we determine the 
matrices representing the exact interpolation functions of effort. The external forces of 'finite element' are supposedly similar 
with the same nature as the internal forces of the same element. 

One of the methods to calculate the external forces of "finite element" is the linearly interpolated from the equilibrium 
equations of the system. Notably in our study efforts are assumed constant at every point of "finite element" and moments 
vary linearly as a function of 

its variables (x and y in case of a plate). Thus, for a triangular plate finite element UK, we obtained the following 
relationships: 



The 



that binds the membrane and bending efforts on any point with the reduced efforts is defined by: 

E m/ } = {N,M} T = {N xx ,N xx ,N xx ,M XX ,M , y M J\ T = [D cmf (^,T])~\{i: r J } (l) 

ie matrix that binds the shear efforts on any point with reduced efforts is defined by: 

{r} = {r„r„} r =[fc„]{s;} e> 
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Where W^j the vector of nodal efforts reduced, D, (<^,T]) and Z? f( are the matrices that represent accurate 

interpolation functions of the efforts membrane bending and shear respectively in the absence of apportionment. The 
stiffness matrix is simply the inverse of the flexibility matrix. 

I ^cmf I an ^ I -Tt ' -^y I are respectively the vector normal forces, effort membrane, bending moments, twisting moment and 
shear forces applied to the cell. 

The direct connection of the finite element provides the stiffness matrix of elementary model in the local 
coordinate expressed by: 

[k<] = [r] t [f; x } 1 [r] (6 ) 

[ p *T =([ f *w)] + [ f *w]Y (7) 

Where F„ x (cmf) and F „ x (cis) are respectively the flexibilities of the matrices membrane combination 
bending and shearing of the plate. I R I is the transition matrix to the system without rigid modes of deformation within five 
degrees of freedom, whose force field is represented by equation (8) and the corresponding displacements <q \ are defined 
(eqt.9): 

[F*"]=[*r{2;} (8) 

{q} = [R]{u} (9) 

With p paq me external force exerted by a plate finite element nodal loads equivalent to the same element and \ U e \ the 
corresponding vector of nodal displacement 1 - and is given by equation (10): 

[u } = {w 0/ , v 0; , w Qi , p Xi , p yi , u 0/ , v 0; , w 0/ , (3 Xi , (3 yj u 0k , v 0k , w Qk , (5 Xk , p yK } (10) 

Remark: In the simple case of a beam with two nodes with three degrees of freedom [23] the force vector corresponds 
exactly to the demands of the nodal finite element beam. 



Flexibility matrices concerning the plates are given by: 

[FfZ{cmf)\ = S m X S H[D cmf ^,rj)'J [H cmf ^ K ,r, K )] X [D^^rj^drj (ii) 

[F^icisaill)] = S IJK "^\\[b ct ] T [H ct {$ K , lK )Y [b ct ]d^d V (12) 
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The matrices H cm A^ K ,T] K ) and \H ct (t^ K ,T] K )\ are matrices named flexibilities membrane bending and shear 

respectively: 

[D^ K ,rj K )]-[H c ^ K ,rj K )]{d cnf ^ K ,rj K )} and fc f H# cf (4^)]R,} 
[H cmf (^ri K )] = r m T H H mf ] ^[HM K ^ K )] = N T h i H n 

with H m =Y j hfl., H f =- £ (z? +1 -Z-)H i and H mf = J /*./;.#. 
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:s\ H, (g K ,T] K ) and I H c (<^ K ,7] K } I respectively represent the stiffness of membrane bending and 
shearing of the cell k of the plate, hi and 27 represent respectively the thickness and position Z layer ;' of the cell, E t and V i 
being respectively the Young's modulus and Poisson's ratio of the corresponding layer, k is the shear correction factor. 
| d-cmf (4k ' 'Hk ) | = | e u : ' e yv ' /^tv ' ^tt ' ^yv ' Xxy \ is me vector of plane deformation, and membrane of curvature 
experienced by a cell, and •! d ct \ = \ y x , y \ is the vector of deformations of the distortion in the planes (x, z) and (y, z). 

3. PRESENTATION OF AN ELEMENT DKT (Discrete Kirchhoff Triangle) 

The DKT element defined in [1J is a finite element with three nodes and three degrees of freedom per node. It is 
considered in this article, as a finite element with three node and five degrees of freedom per node. 

The rotations /3 x , B are interpolated in a parabolic manner and the transverse displacements U , V Q , W Q are interpolated 
in a linear manner [1,2]: 

p x = ^N ( j3 Xi + J P x a k ,p y =Y,N i j3 yj + J P y a k , P Xt = P k C k and P v , = P k S k 
u = y\N t u , v = V N t v , w = V N t w 



Where C k , S k are the direction cosines, k is the middle of respective sides of the triangle, and are given by the side 

if. C k = (Xj -x t )IL k , S k = {y j - y t )IL k and L k = J{x j -xf +(y. - yf 

Where n is the number of nodes of the finite element, in the case of a triangular element Tl — 3 and functions N t and P k 
are given by [1,2]: 

N l= A = l-%-r/, N 2 =%, N 3 =r/, P 4 = 4%A, P 5 =4^rj and P 6 = ArjX 

The expression of CC k according to the nodal variables of nodes i and / is [1]: 
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- ( W . - Wj ) — (C k j3 Xi + S k j3 yj + C k /3 Xj + S k j3 yj ) 
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III. ANALYSIS OF A UNIFORM PLATE WITH DKT AND FLEXIBILITY (FLX) 

At first glance, the figure 2 represents the results obtained with FLX as we analyze a homogeneous square plate 
subjected to uniform load simply supported or built on the contour, for dilTcienf slenderness L/h (5 to 1000). The plate is 
meshed with 128 (N = 8) rectangular isosceles elements (§ 2.2). The results are virtually identical with those obtained with 
DST and Q4y [1] for the recessed plate (Figure 2. a). For the simply supported plate there appeared an error of 
about -0.5% (Figure.2.b). 
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b) Simply supported plate 
Figure 2-homogeneous square plate with load uniform. Bending in the centre based on L/h 

(D = Eh 3 /I2(l-v 2 ) ,v = 0.3 ,k =5/6) 
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In a second step, we describe in Figures 3 to 7 some results [1,2] and on the analysis of a homogeneous square plate 
subjected to a concentrated load at the center or simply supported and built on the contour. A quarter of the plate is meshed 
with 2, 8, 32, 128 (N = 1, 2, 4, 8) rectangular isosceles DKT elements (§ 2.3) and FLX (§ 2.2). These elements have five 
degrees of freedom per node and are of Kirchhoff ( no iranw erse shear energ) , the results rue independent of UK) for DKT 
elements and flexible elements taking into account of transverse shear (FLX) for LI h — 24 (figures 6 and 7). Figures 4 
and 5 we presents the results obtained with FLX as we analyze a homogeneous square plate subjected to concentrated load 
simply supported or built on the contour, for different slenderness L/h (5 to 1000) for both types of mesh (there is a thin 
outlook of influence mesh ). The plate is meshed with 2, 8, 32 and 128 (N = 1, 2, 4,8) rectangular isosceles FLX elements 
(§ 2.2). Convergence can be seen for N = 8, that is to say, for a mesh of 128 elements, we observe a occurrence of an error, 
for the clamped plate, in the order of —0.5% mesh A (Figure 4. a) and —0.56% mesh B (Figure 5. a) and simply 
supported plate —0.3% mesh A (Figure 4.b) and —0.31% mesh B (Figure 5.b). In Figure 6, we provide the percentage 
error of the deflection at the center depending upon 'N' number of divisions per half side. There is 

with FLX (FLX model is a consistent shift, the total potential energy E p > E p and as E p ■ 

that(w c ) <(w c ) 
convergence of DKT ( 
elements of the type DKT and FLX. The convergence of the moment \M x in the 

reaction concentrated in the corner (= Y1M ) in case of simply supported plate are presented in Figure 7 for both types 

of meshes and for DKT and FLX. (Calculations of efforts have been made directly to the nodes peaks followed by an 
average if the node is shared by two elements). There is a fairlj rapid convergence, an influence of models and an orientation 
of the mesh. 
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• Symmetry conditions: 
J3 X =0 onCA;onCD^=0 

• Boundary conditions: 

- Recess ;w = p x =/3 y = on ABD 

- Support simple: w = f3 x = on AB, 
w = a = on BD 



Meshes considered: N= 
Case N = 2 




mesh A 

Kirchhoff solution for a concentrated load P: 



(D = Eh 3 /12(l-v 2 );v = 0.3) 

Recess: W =5.6 xlO" 3 PL 2 / D and \M r | = 0.1257P 
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Simple Support: W c = 1 1.6 X 10" 3 PL 2 I D and /? = J2M vy | = 0.1219P 
Figure 3-square plate under concentrated load. Data 
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simple Support 

Figure 5-square plates with concentrated load at center built and simply supported. Error for DKT and FLX 
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simple Support 
Figure 6-square plates with concentrated load at center built and simply supported. Error Ot 
the corner for DKT and FLX 

IV. CONCLUSION 

The flexibility method developed with a linear interpolation (interpolation functions of the first order) and of way 
independently of the transverse displacements and rotations, solves the problem related to the phenomenon by 
shear locking. In the case of multicellular multilayer finite element, we observe that the method of flexibility, 
which is a model monotone convergence, converges quickly enough for a plate structure. In this paper we have 
presented the results for the analysis of a square plate in bending under load concentrated at the center, simply 
supported on the contour or clamped while highlighting the influence of the mesh on different slenderness L / h 
(Figures 4 and 5: arrow report W c I W k ). We also presented results on an analysis of a square plate subjected to a 
uniform load, clamped or simply supported on the contour (Figure 2 ). The percentage error appeared in Figures 4, 
5, 6 and 7 and that can be translated by the phenomenon of blocking is reduced (becomes negligible) by increasing 
the number of elements this allows us to confirm the reliability of the method on solving the problem of shear 
locking. In the following work (in a future article) we present the results at predictive calculation of the 
performance of bearing subject to the sails seismic behavior by numerical simulation coupled with a damage 
model by comparison with experimental results and by adopting a damage model for multicellular multilayer 
finite element . 
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